Atomic-scale surface demixing in a eutectic liquid BiSn alloy 
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Resonant x-ray reflectivity of the surface of the liquid phase of the Bi43Sn57 eutectic alloy reveals 
atomic-scale demixing extending over three near-surface atomic layers. Due to the absence of un- 
derlying atomic lattice which typically defines adsorption in crystalline alloys, studies of adsorption 
in liquid alloys provide unique insight on interatomic interactions at the surface. The observed com- 
position modulation could be accounted for quantitatively by the Defay-Prigogine and Strohl-King 
multilayer extensions of the single-layer Gibbs model, revealing a near-surface domination of the 
attractive Bi-Sn interaction over the entropy. 

PACS numbers: 61.20.-p, 61.10.-i, 68.03.-g 
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The widely-accepted Gibbs adsorption rule pre- 
dicts the surface segregation of the lower surface en- 
ergy component of a binary mixture. Liquid metals are 
ideal systems for studying Gibbs adsorption due to the 
nearly spherical shape of interacting particles, relative 
simplicity of the short-range interactions and the avail- 
ability of bulk thermodynamic data for many binary al- 
loys. While certain aspects of Gibbs theory can be tested 
through macroscopic measurements of surface tension or 
adsorption isotherms, very few direct measurements of 
the atomic-scale composition profiles of the liquid-vapor 
interface were reported 0,0,3- In addition to funda- 
mental questions related to surface thermodynamics of 
binary liquids, BiSn-based alloys have been widely stud- 
ied as substitutes for Pb-based low-melting solders 
Thus, understanding their wetting, spreading, alloying, 
reactivity and other surface-related properties is of great 
practical importance. Moreover, interfacial phenomena 
dominate the properties of the increasingly important 
class of nanoscale materials, as demonstrated recently in 
studies of the liquid-solid phase stability of nanometer- 
sized BiSn particles 

Synchrotron-based x-ray reflectivity (XR) can mea- 
sure the surface-normal density profile of a liquid with 
Angstrom-scale resolution. Over the last decade XR re- 
vealed the long-predicted surface-induced atomic layer- 
ing at the liquid- vapor interface for a number of elemen- 
tal liquid metals Q, S S E El. Resonant XR near 
an absorption edge resolved the density profile of each 
component in Gain 0], HgAu and Biln liquid bi- 
nary alloys. The enhancement of the concentration of the 
low-surface-tension component was invariably found to 
be confined to the topmost surface monolayer, with sub- 
sequent layers having the composition of the bulk, in ac- 
cord with the simplest, and widely used, interpretation of 
the Gibbs rule. By contrast, we find here an atomic-scale 
phase separation extending over at least three atomic 



layers. This is unexpected, considering the near-perfect- 
solution nature of the Bi43Sn57 alloy in the bulk p^ll3|. 

A hquid Bi43Sn57 sample (99.99% purity, Alfa Aesar) 
was prepared under UHV conditions (P< 10~^ Torr). 
Atoniically clean liquid surfaces were obtained by me- 
chanical scraping and Ar"*" ion sputtering, as described 
previously 0, Ej, 0| . Measurements were done using 
the liquid surface diffractometer at the ChemMatCARS 
beamlinc. Advanced Photon Source, Argonne National 
Laboratory at a sample temperature of T = 142 °C, 
4 °C above the Bi43Sn57 alloy's eutectic temperature, 
Te = 138 °C. 

The reflected intensity fraction, R{qz), of an x-ray 
beam impinging on a structured liquid surface at a graz- 
ing angle a, is given by the Born approximation as: 



R{q.) = Rpiqz) ■ mq^)\^ ■ CWiq,) 



(1) 



where — (47r/A) sin a, A is the x-ray wavelength, 
Rpiqz) is the Fresnel XR of an ideally abrupt and flat 
surface, CW{qz) is due to thermal surface ca pilla ry waves 
0,0], and the surface's structure factor is |lfil|: 
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Here z is the surface-normal axis, poo and p{z) are the 
bulk and surface electron densities, respectively, and (..) 
denotes surface-parallel averaging. As Rf is a univer- 
sal function depending only on the known critical angle 
for total external reflection, and CW{qz) is known accu- 
rately from capillary wave theory, the intrinsic density 
profile, {p{z)), is obtained by computer fitting the mea- 
sured Rigz) by a physically motivated model described 
below [lOj. 

The (forward) atomic scattering factor of a Z-electron 
atom varies with energy as^: Z' = Z + f {E)-if" (E), 
where f'{E) and f"{E) = p[E)\/ {Att), are the real and 
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FIG. 1: Dispersion corrections (a) f'{E) and (b) f"{E) of 
Bi near the L3 absorption edge at El3 = 13.418 keV. The 
right scale of (b) is the electron density contrast AZ' / Z'^ = 
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FIG. 2: XR measured at the indicated energies, with fits by 
the three-layer model (lines). The dashed line is the XR of a 
uniform-composition surface. Inset: The E=12.00 keV mea- 
sured R/Rf with fits by the three models discussed in the 
text (lines). Error bars are smaller than the symbols' size. 



imaginary parts of the dispersion correction, and fJ-{E) is 
the photoelectric absorption coefficient. 

The effect of /" on the analysis can be neglected and 
the changes in /' are significant only near an absorp- 
tion edge. Fig. ^a) shows f"{E) near the Bi L3 edge 
as obtained from an absorption measurement in a Bi 
foil. Fig. ^b) is the corresponding f'{E), calculated 
from f"(E) using the Kramers-Kronig relation 0| . Both 
agree well with theory [isl Il9| . The composition depen- 
dence of {p{z)) was obtained by fitting the measured XR 
by the distorted crystal (DC) model for a layered liquid 
surface 0,|i|: 

{p{z)) ^ e-(--"'^)V-^ / Z'^, - Z', .c- m 

Poo ~ V2^ajd I' + Z'^ ) c 

The progressive increase in the Gaussian width parame- 
ter cr^ = (Tg -I- (n — 1)ct^ with increasing layer number n 
describes the layering amplitude's decay below the sur- 
face 0. The layer spacing d is kept constant in this 
model due to similarity in size between Bi and Sn atoms. 
The bulk's average effective electron number per atom 
is Z'^ = xZg^ + (1 — x)Zg^, and (5„ = x'^ ~ x is dif- 
ference in the Bi fraction between the rt-th layer, x[^, 
and the bulk, x. The corresponding atomic densities, 
c„ and c, are determined from the atomic volumes VBi 
and vsn- Cnx'^vsi + c„(l — x'„)vsn = 1- The contrast, 
{Z'g^ — Zg^)/Z'^, varies strongly near the edge due to 
the variation of Z'^-: from 0.43 at E=12.00 keV to 0.27 
at E=13.418 keV (right axis in Fig.[Tl). This is the basis 
for the resonant XR method which allows to separate out 
the density profiles of the two species 0, . 

Fig. 121 shows Fresnel-normalized XRs R(qz)/ Rpiqz) 



measured near the Bi L3 edge at the four energies marked 
by triangles in Fig. ^b). The dashed line is calculated 
from the DC model for a layered interface assuming a 
uniform composition ((5„ = 0). The strong enhancement 
of the measured R{qz)/ Rf{<1z) over this line, evidenced 
by the peak at ~ 1.0 A~^, and the strong energy- 
dependence of the low-qz reflectivity, clearly indicate a 
significant surface segregation of Bi, and its variation 
with z. 

Three fits of the data by the DC model, Eq. O were 
carried out, assuming that only one {Si ^ 0, Sn>2 = 0), 
two (Ji,2 7^ 0, 5n>3 = 0), or three {61^2.3 7^ 0, i5„>4 = 0) 
surface layers deviate from the bulk composition. All fits 
employed d = 2.90 A, ctq = 0.30 A and = 0.57 A, de- 
rived from the energy-independent position, shape and 
intensity of the layering peak at q^ = 2.0 A^^. The mea- 
sured R{qz) / Rf{<Iz) of all four energies were fitted simul- 
taneously, using the experimentally determined f'{E). 

Fig.|21exhibits an excellent agreement of the three-layer 
model (solid lines) with the measurements, but a very 
poor agreement for the one- and two-layer models (inset) . 
Table m lists the best-fit values of x'^ and (5^** = x',^ ~ x 
and the corresponding 95% non-linear confidence inter- 
vals Y{x'„) and F((5^**) determined from a six-parameter 
support plane analysis jl^. The most striking result is 
the non-monotonic deviation Sn of Bi from the 43% bulk 
value, showing an enhanced composition of 96% and 53% 
in the first and third layers, and depletion down to 25% 
in the second layer (see Fig. 3). Beyond the third layer 
entropy effects dominate the Gibbs adsorption and the 
layer and bulk concentrations can not be distinguished. 
While demixing has not been previously reported in liq- 
uid alloys, similar decaying oscillatory composition pro- 
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FIG. 3: Electron density profiles as derived from the fits to 
the reflectivities shown in Fig. |2| Inset: the bulk-normalized 
difi'erences in electron density of Bi and Sn, (psi — psn)/ poo- 



files were discovered in several crystalline alloys such as 
CU3AU [ilj . However, the properties, formation mech- 
anism, and strong temperature dependence of the com- 
position modulations in CuaAu alloys were found to be 
intimately related to, and largely dominated by, the long- 
range fee order in the bulk crystal, and the severe packing 
strains resulting from the big mismatch in the atomic di- 
ameter of the two components. As none of these exist 
in our liquid alloy, surface-induced segregation and the 
Gibbs rule can be studied in a pure short-range-order in- 
teraction context, free from the complicating influence, 
or even dominance, of other effects. We now compare 
our experimental observations with theory. 

Guggenheim's application of Gibbs theory Q to 
regular solutions assumes the surface segregation to be 
restricted to a single surface monolayer. Assuming p 
nearest neighbors for each bulk atom in a layered lat- 
tice model, Ip are within, and mp are in the adjacent, 
layers. For a close-packed lattice, for example, p = 12, 
I = 0.5 and m = 0.25. The surface tension of the regular 
solution, 7_AB, follows from those of the pure components, 
7A and 75, as jH]: 



lAB = IB 

kT 
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Here, x and (1 — x) are the bulk concentrations of atoms 
A (Bi) and B (Sn), while x' ^ x and (1 — x') are the corre- 
sponding surface concentrations, a a and as are the two 
atomic areas, and to = 2ujab — ^aa — ^bb is the interac- 
tion parameter, defined by the A-B, A-A and B-B atomic 
interaction energies. Extrapolated down to T = 142 °C, 
7Bi = 398 mN/m and 7s„ = 567 mN/m, while asi and 
asn are calculated from the atomic radii rsi = 1.70 A 
and rsn = 1.62 A assuming hexagonal close packing [23| . 
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TABLE I: Density model parameters x'^ and 5^'' — x'„ — x, 
and confidence intervals Y(x'„) and Y{5^^^) obtained from the 
three-layer model fits compared to theoretical and Sf^^ 
derived from the Defay-Prigogine and Strohl-King models. 



Treating Bi43Sn57 as a perfect solution [uj/kT = 0), the 
Gibbs theorem, Eq. ^ yields ^ab = 444 mN/m and 
x' = 0.904, below the experimental value a;i(Bi) in Table 
m However, assuming a regular solution behavior with 
uj/kT ^ 1 yields jab = 432 mN/m, and x' = 0.941, 
which agrees very well with the experimentally derived 
^^(Bi) in Table Both 'jab agree well with experiment 
and theory j^l ■ Note that jab and x' are only weakly de- 
pendent on w/ kT due primarily to the logarithmic func- 
tional behavior and large surface tension difference of the 
two components, asnJSn — o-BilBi ~ 2 kT. This intro- 
duces a large uncertainty of co/kT calculated from mea- 
surements of surface tension or surface monolayer compo- 
sition. Resonant XR measurements of sub-surface layer 
composition therefore present a unique opportunity to 
probe the nature of atomic interactions at the surface. 

In spite of the good agreement above, confining the 
surface excess to a single monolayer is correct for perfect 
solutions only, but not for our case of a regular solution, 
as Defay and Prigogine [2^ point out. They provide a 
correction for regular solutions, where the surface excess 
extends over two layers, the 7^5 values above do not 
change significantly and the layers' (5„ are related by: 
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1 + S2/X 



2uj , 
kT ^ 



2ijjm 



1 - (52/(1 -x) kT"^ kT 
Expanding Eq. [S] to first order in 82 '■ 
2ujmx{l ~ x)5i 



(<5i - 282) = 0. (5) 



52^ 



kT -2ujlx{l 



(6) 



For nearly perfect solutions [uj/kT ^ 1) Eq. [S] yields 
a negligible 62'- Q < 82 ^ 5i. For ujjkT > 1, how- 
ever, 82 and Si are of opposite signs and \S2\ may be- 
come comparable to |i5i|. This prediction is qualitatively 
consistent with the demixing observed here. For exam- 
ple, when Lo/kT 3> 1, Eq. El can be simplified further: 

52 ~ —{m/l)5i. For Bi43Sn57, m/l w 0.5 and the Gibbs- 
predicted x[ = 0.90 (or Si = 0.47) yields 62 -0.23, 

53 = 0.12 and S4 ~ —0.06 |2^. These values, shown as 
S^^ in Tabled agree well with (5^** obtained from the 
three-layer model fits. The smallest value of the inter- 
action parameter oj /kT for which satisfactory agreement 
with the Defay-Prigogine model could be obtained (by 
treating m as an adjustable parameter) is ujjkT = 2.3. 
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Strohl and King [23| suggest a multilayer, multicompo- 
nent model, where no expansion is used, and x'^ are ob- 
tained iteratively, until convergence to a self-consistent 
composition profile is reached. A good agreement of this 
theory with our BiSn measurements is obtained when 
Lo/kT = f.O - 1.7. Typical 5^^ values are listed in 
Table HJ As observed, the Strohl-King model provides 
composition profiles very similar to those of the Defay- 
Prigogine model, albeit with slightly different 6n values, 
thus supporting our overall conclusions. 

Theoretically, lj and the enthalpy of mixing, AiJ,„, 
are related by = l^Hm/[x{l — x)]. In practice, how- 
ever, bulk thermodynamic quantities were often found to 
yield inaccurate values for surface quantities. For exam- 
ple, organic jl^ and metallic [l^ mixtures exhibit sig- 
nificant disagreements between uj values derived empiri- 
cally from surface tension measurements and from bulk 
calorimetry. For BiSn, reported values of Affm range 
from endothermic values of 80 to 140 J/mol [13 to an 
exothermic value of -180 J/mol These values lead 

to \ijo/kT\ < 0.2, i.e. an almost perfect solution, and an 
insignificant |^2| < 0.01. On the other hand, the value 
of w/fcT w 10 that we previously found necessary to ac- 
count for the observed 35% Bi concentration enhance- 
ment at the surface monolayer at the Biln eutectic is of 
the same order of magnitude as the value we find nec- 
essary to account for the present observation of surface 
segregation in BiSn, cu/kT ^ 1.0 — 2.3. Unfortunately 
we do not have an explanation for the origin of the dis- 
crepancy in the values of uj/kT and this suggests an ur- 
gent need for both further theoretical studies of surface 
demixing as well as experimental investigations of simi- 
lar effects in other binary alloys. In particular, the BiSn 
system appears to be the only liquid alloy for which clear 
evidence for multilayer surface demixing has been found. 
The case for new studies is strongly reinforced by the 
existence of a growing class of surface-induced ordering 
phenomena that have been observed in metallic liquids. 
In addition to the surface demixing reported here, these 
include laye ring relaxation |0, segrega- 

tion 0, l^iM I13| , wetting transitions 0, |^ , and surface 
freezing j32l | . Finally, there is a basic unresolved question 
of whether the surfaces of liquid metals are fundamen- 
tally different from those of non- metallic liquids [ssj . 
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